TUM-HEP-667/07 
MPP-2007-43 



Correlations between e'/e and Rare K Decays 
in the Littlest Higgs Model with T-Parity 

Monika Blanke ab , Andrzej J. Buras a , Stefan Recksiegel a , 
Cecilia Tarantino a and Selma Uhlig a 

a Physik Department, Technische Universitat Miinchen, D-85748 Garching, Germany 
b Max-Planck-Institut fur Physik (Werner-Heisenberg-Institut) , 
D-80805 Miinchen, Germany 



Abstract 

We calculate the CP-violating ratio e'/e in the Littlest Higgs model with T-parity 
(LHT) and investigate its correlations with the branching ratios for Kl — > k^vv, 
Kl — > ir°£ + £~ and K + — > it^vv. The resulting correlations are rather strong in the 
case of Kl decays, but less pronounced in the case of K + — > Tt + vv. Unfortunately, 
they are subject to large hadronic uncertainties present in e'/e, whose theoretical 
prediction in the Standard Model (SM) is reviewed and updated here. With the 
matrix elements of Qq (gluon penguin) and Qs (electroweak penguin) evaluated in 
the large- TV limit and m^ s (2GeV) = 100 MeV from lattice QCD, (e'/e) S M turns 
out to be close to the data so that significant departures of Bt{Kl — > tt^vv) and 
Br(Ki — > tt°£ + £~) from the SM expectations are unlikely, while Br(K + — » it + ui>) 
can be enhanced even by a factor 5. On the other hand, modest departures of the 
relevant hadronic matrix elements from their large- N values allow for a consistent 
description of e'/e within the LHT model accompanied by large enhancements 
of Br{Ki — > tfivv) and Br{Ki — > tt £ + £~), but only modest enhancements of 
Br(K + -> tt + vv). 



Note added 



An additional contribution to the Z penguin in the Littlest Higgs model with T-parity 
has been pointed out in [1,2], which has been overlooked in the present analysis. This 
contribution leads to the cancellation of the left-over quadratic divergence in the cal- 
culation of some rare decay amplitudes. Instead of presenting separate errata to the 
present work and our papers [3-6] partially affected by this omission, we have presented 
a corrected and updated analysis of flavour changing neutral current processes in the 
Littlest Higgs model with T-parity in [7]. 



1 Introduction 



Flavour Changing Neutral Current (FCNC) processes provide a powerful tool for testing 
the Standard Model (SM) and its extensions. Of particular interest are the four rare kaon 
decays K L -» 7r°z/z/, K + — > tt^uu, K l — > ir°e + e~ and K L — > n°Lt + Li~. Their branching 
ratios are strongly suppressed within the SM and consequently can be largely modified 
by New Physics (NP) contributions. 

Extensive analyses of these decays in the MSSM [8], the Littlest Higgs model with 
T-parity (LHT) [3], general models with enhanced Z-penguin contributions [9] and Z'- 
models [10] have shown that in the presence of new sources of flavour and CP-violation 
beyond those present in the Minimal Flavour Violation (MFV) framework [11,12], en- 
hancements of Bt(Kl — > Ti Q vv) by an order of magnitude and of the other branching 
ratios by up to a factor 5 are still possible. 

On the other hand, as pointed out in [13] and analyzed in more detail within the 
MSSM in [14], the enhancements of the rare decay branching ratios in question could be 
bounded in principle by the value of e'/e that measures the ratio of the direct and indirect 
CP-violating contributions to Kl — > tttt. The reason is very simple. The electroweak 
penguin and box diagrams that enter the evaluation of the rare decay branching ratios 
in question have also considerable impact on the ratio e'/e so that, in a given model, 
specific correlations between e'/e and the branching ratios for rare K decays exist. 

Unfortunately, whereas the branching ratios of K — > nvv decays are theoretically 
very clean [15] and those of K L — > tt°£ + £~ are subject to only moderate theoretical 
uncertainties [16], this is not the case for the ratio e'/e that is affected by large hadronic 
uncertainties. 

Indeed, whereas the Wilson coefficients of the local operators entering the evaluation 
of e'/e are known [17-22] at the NLO level in QCD and QED renormalization group 
improved perturbation theory, the hadronic matrix elements of these operators are still 
only poorly known. 1 Therefore the predictions for e'/e in the SM and its extensions have 
very large theoretical uncertainties. 

In spite of this unsatisfactory situation and in view of future improvements in the 
evaluation of the relevant hadronic matrix elements by lattice QCD or large- N methods, 
we think that it is important to analyze the correlations between e'/e and rare kaon 
decays in specific extensions of the SM, where large enhancements of the rare decay 
branching ratios have been found. Certainly, the result of such an exercise will sensitively 
depend on the values of the hadronic parameters present in e'/e, but the mere fact that 
such correlations exist will hopefully motivate further non-perturbative studies. 

The main goal of the present paper is the calculation of e'/e within the LHT model 

1 Latest reviews can be found in [23-27]. 
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[28,29] and the investigation of its correlations with the four rare kaon decays in question, 
for a given set of the non-perturbative parameters entering e'/e. To this end we will 
apply a useful parameterization of e'/e proposed in [23] that automatically takes into 
account all renormalization group effects from scales below m t and expresses the hadronic 
uncertainties in terms of the two parameters P 6 and R 8 corresponding to the dominant 
QCD and electroweak penguin operators, respectively. 

In [3] very sharp and theoretically clean correlations between the decays K L — > ir°v9, 
K + — > ix+vv and K L — > have been found in the LHT model, subject mainly to 

a discrete ambiguity present in the correlation between Kl — » n®vv and K + — > it + vv. 
It is therefore sufficient to establish the correlations between e'/e and K L — > -k^vv and 
between e'/e and K + — > n+uu in order to get an idea about all correlations. 

Our paper is organized as follows. In Section 2 we briefly review the status of e'/e 
in the SM, investigate the relevant theoretical and parametric uncertainties and provide 
a numerical update of [23]. In Section 3 we present the basic formulae for e'/e in a 
generic model with new complex phases but no new operators relative to the SM, in 
terms of the short distance functions X, Y, Z and E that contain both SM and NP 
contributions. It turns out that in the LHT model the functions X, Y and Z can directly 
be obtained from our previous analysis [3]. The function E that plays a subdominant 
role in e'/e, is calculated for completeness here for the first time in the LHT model. 
In Section 4 we evaluate e'/e in the LHT model scanning over its parameters and for 
various values of Rq and R$. The main results of this section are the correlations between 
e'/e and Bt(Kl — > tx q uu) and between e'/e and Br(K + — > -k+vu), that illustrate the 
very important role of e'/e in bounding the enhancements of rare K decay branching 
ratios provided the non-perturbative parameters Rq and i?§ are accurately known. We 
conclude in Section 5. 



2 e'/e in the SM 
2.1 Basic Formula 

Before analyzing e'/e within the LHT model, it will be instructive to have a brief look 
at this ratio within the SM and investigate the relevant theoretical and parametric un- 
certainties that have to be taken into account also in the case of the LHT model. This 
will also allow us to update the analysis of [23]. 

The formula for e'/e of [23] is given in the SM as follows: 



e' 



Im(At) • [P + P E E (x t ) + PxX (x t ) + P Y Y {x t ) + P z Z (x t )j (2.1) 



2 



with At = V t * s V td and Xt = mf/M^. The short distance physics is described by the 
loop functions E (x t ), X (x t ), Y (x t ) and Z (x t ), for which explicit expressions can be 
found in [30]. On the other hand, the Pi encode information about the physics at scales 
ji < 0(m t , M\y), and are given in terms of the hadronic parameters 



as follows: 



121 MeV 



m s (m c ) + m d (m c 



D _ R (3/2) 



rf + rf ^ + rfR 8 . (2.3) 



121 MeV 



m s (m c ) + m d (m c 



(2.2) 



The coefficients r\°\ rf^ and rf^ enclose information on the Wilson-coefficient functions 
of the AS* = 1 weak effective Hamiltonian at the next-to-leading order [30]. Their nu- 
merical values for different choices of at a — m c in the NDR renormalization scheme 

MS r° c 

can be found in [23]. The numerical values of the Pi are sensitive functions of i? 6 and R s , 
as well as of or equivalently a s (M z ). The values = 310, 340, 370 MeV consid- 
ered in [23] and by us correspond to the three-loop values a s (M z ) = 0.119,0.121,0.123, 
respectively. The two- loop formula for the strong coupling constant, instead, provides 
a s (M z ) = 0.117,0.119,0.121. Although three-loop values are quoted by the PDG [31,32], 
we use in the present analysis the two-loop values, as the Wilson coefficients entering 
e'/e are known at the NLO only. 

2.2 Status of B^ /2) and B ( 8 3/2) from Lattice QCD 

The hadronic parameters B^ 2 ^ and B^ 2 ^ represent the matrix elements of the domi- 
nant QCD penguin operator Qq and the dominant EW penguin operator Q§, respectively. 
They are the main source of uncertainty in the determination of e'/e and, hence, cal- 
culating (7T7r| Qq$\K) reliably represents a theoretical challenge for the non-perturbative 
methods like lattice QCD and large- N. The large- approach will be referred to below 
while we focus here on the status of lattice studies relevant for e'/e. The lattice calcu- 
lation of the Qq matrix element is particularly delicate. Golterman and Pallante [33], 
indeed, have pointed out that there is a serious ambiguity in the lattice version of left- 
right QCD penguin operators, like Q 6 , because the flavour group in (partially) quenched 
QCD is not SU(3) but SU(3 + Nf\3) where N f is the number of sea quark flavours. 
It turns out that the ambiguity in Q 6 has such a large effect on e'/e that it can even 
flip its sign in quenched QCD [34-36]. Moreover, the same problem affects the left-left 
QCD penguin operator Q 4 with a sub- leading effect in e'/e [37]. On the other hand, 
the lattice calculation of the Q 8 matrix element is more reliable, although challenging 
as well and still affected by an uncertainty of 10 20%. Two independent approaches 
have been used. In the indirect approach, one calculates the hadronic matrix elements 
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of K — > 7r and K — > and reconstructs K — > txtx amplitudes using chiral perturba- 
tion theory. This method, relatively easy and computationally cheap, has been widely 
used [34,38,39], but it only works in leading order chiral perturbation theory. In the 
direct approach, instead, one calculates directly the K — > txtx matrix elements with the 
final state pions carrying a physical momentum. The difficulty of this method is repre- 
sented by the Maiani- Testa "no-go theorem" [40]: one can not obtain K — > Tx(p)Tx(—p) 
but only K — > 7r(0)7r(0) on the lattice, where |7r(0)7r(0)) is the ground state of two pions 
with periodic boundary condition in the spatial direction. Various methods have been 
proposed to get around the Maiani- Testa theorem. Luscher and Wolff [41] proposed a 
diagonalization method, based on a computationally expensive calculation of correlators 
with non-zero pion momentum. Another possibility consists in modifying the boundary 
condition for the pions [42], thus providing a finite momentum to the ground state of 
tx ± . A different approach was elaborated by Lellouch and Luscher [43], based on an 
excited state fit to extract the \7x(p)7x(— p)) state that appears in a finite volume where 
the spectrum of two-particle states is discrete, and on a formula for connecting the decay 
measured in a finite volume to the infinite volume result, in the center of mass (CM) 
frame. This technique, however, is challenging due to the need to extract the excited 
state. An alternative and promising method is to work with a kinematic setup for which 
the final state of interest is also the lowest energy state. This has been done in [44-46] 
by working in the moving (LAB) frame, i.e. calculating {7x(P)7x(0)\Qs\K (P)) and then 
converting the result from the finite to the infinite volume, using the Lellouch-Luscher 
formula [43]. An important theoretical advance of the last years is the derivation of a 
relationship similar to the Lellouch-Luscher formula but valid in the LAB frame [47, 48] 
that may improve the accuracy of the LAB-frame method, as shown in a preliminary 
calculation with domain wall fermions [49]. 

2.3 Comparison between SM Prediction and Experimental Data 

On the experimental side, the world average based on the latest results from NA48 [50] 
and KTeV [51] and previous results from NA31 [52] and E731 [53] reads 

e'/e= (16.7 ± 1.6) • 10" 4 . (2.4) 

While several analyses made in recent years within the SM found results that are com- 
patible with (2.4), it is fair to say that the large hadronic uncertainties in the coefficients 
Pi still allow for sizeable NP contributions. The relevant list of references can be found 
in [23-27]. 

In [23] an agreement of the SM with (2.4) has been found for (i? 6 , R$) = (1.2, 1.0) and 
A^ = (340 ± 30) MeV. Meanwhile the value of m t decreased and the value of Im(A t ) 
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increased. Consequently for R§ = R$ = 1.0, corresponding to the large-iV approach 
of [54] with B ( l /2) = Bf /2) = 1, and mf(2 GeV) = 100 MeV from lattice QCD [31, 55], 2 
acceptable agreement with (2.4) can be obtained, provided A^- > 340 MeV. Indeed in 
this case we find for Agi = 340 MeV 

P = 15.962, P x = 0.597, P Y = 0.519 , P z = -12.416, P E = -1.226, (2.5) 

and choosing lm(\ t ) = 1.38 • 10~ 4 , obtained by the UTfit collaboration [57], the result 

(eVe)sM = 12.3 • 10" 4 (2.6) 

which is a bit lower than the value in (2.4). For A^ = 370 MeV we find, on the other 
hand, 

(e7e) S M = 13.5-10- 4 , (2.7) 

within 2(7 from the central value in (2.4). A slight decrease of the m^ s (2GeV) value 
(see Table 1) would result in an improved agreement with the data. 

We would like to emphasize, then, that with Im(A t ) = 1.69 ■ 10~ 4 , obtained from the 
tree level determination of the CKM parameters, the values in (2.6) and (2.7) increase 
to 

(£7e) SM = 15.3-10- 4 , (2.8) 

and 

(e'/e) S M = 16.7 -10" 4 , (2.9) 

so that even for A^ = 340 MeV and mf E (2 GeV) = 100 MeV a good agreement with 
the data can be obtained. 

As a preparation for the analysis of e'/e in the LHT model we show in Fig. 1 the 
values of (e'/e) SM for three different choices of (R 6 ,R 8 ) = (1.0, 1.0), (1.5,0.8), (2.0, 1.0), 
different values of A^L and the two values for lm(X t ) considered above. 

The main messages from Fig. 1 and (2.6)-(2.9) are: 

• (e'/e)sM has a visible dependence on the values chosen for lm(X t ) and for A^-, but 
these dependences amount only to about 10 -j- 20%, which is comparable to the 
experimental error in (2.4). 

• (e'/e)sM depends very strongly on the values of R e and i? 8 , and the choices (1.5, 0.8) 
and (2.0, 1.0) give values for (e'/e) SM that clearly are in disagreement with the 
data for the full range of A^ and lm(X t ) considered by us. For instance for 
A^ = 340 MeV and the UTfit value of Im(A t ) one finds (e'/e) S M = 26.3 • 10~ 4 
and (e'/e) S M = 36.5- 10~ 4 for (R 6 ,R 8 ) = (1.5,0.8) and (R 6 ,R 8 ) = (2.0,1.0), 
respectively. 

2 Similar results are found from QCD sum rules [56]. 
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Figure 1: ( y e'/e) SM for three different choices of (R 6 , R 8 ) = (1.0, 1.0), (1.5, 0.8), (2.0, 1.0) 
and different values of A™ = 310, 340, 370 MeV. The values obtained with the UTfit 
value for lm(X t ) UT -f it = 1.38 • 10~ 4 are marked with red diamonds, while those with the 
tree level value lm(X t ) tree = 1.69 ■ 10~ 4 are marked with blue stars. The shaded area 
represents the experimental result in (2.4). 

• Significant although smaller departures of (Rq, R&) from (1.0, 1.0) and therefore 
of e'/e from the data could also occur for B^ 2 ^ = B^ 2 ^ = 1, as obtained from 
the large- N approach of [54] , and values of the strange quark mass deviating from 
mf (2 GeV) = 100 MeV by the present 10 20% lattice uncertainty (see Table 1). 



mf s (2GeV) 


80 MeV 


90 MeV 


100 MeV 


110 MeV 


120 MeV 


(R 6 , R 8 ) 


(1.5,1.5) 


(1.2,1.2) 


(1.0,1.0) 


(0.8,0.8) 


(0.7,0.7) 



Table 1: Choices for the strange quark mass within present lattice uncertainties and 
corresponding values for the hadronic parameters (Rq, R$). The small down quark mass 
has a minor impact and its value is fixed to m^ s (2GeV) = 5 MeV [31]. The variation 
of entering the quark mass running represents a small effect as well and its value is 
fixed to A^ = 340 MeV. 

As reviewed in [23], Rg = 1.0 ± 0.2 is obtained in various approaches. Unfortunately 
the value of Rq is very uncertain. For instance in the large- N approach of [58, 59] 
values for i? 6 significantly higher than 1 have been found. In particular [58] reports 
Rq = 2.2 ± 0.4 and i?§ = 1.1 ± 0.3. On the other hand, while the lattice values of i?§ are 
compatible with 1 [38,44], they are lower than unity for i? 6 [34,35]. 
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3 e'/e in the LHT Model 



The LHT model [28,29] belongs to the class of Little Higgs models [60], where the 
little hierarchy problem is solved by a naturally light Higgs, identified with a Nambu- 
Goldstone boson of a spontaneously broken global symmetry. In the LHT model the 
global group S77(5) is spontaneously broken into 5*0(5) at the scale / « 0(1 TeV) 
and the electroweak sector of the SM is embedded in an SU(5)/SO(5) non-linear sigma 
model. Gauge and Yukawa Higgs interactions are introduced by gauging the subgroup 
of £77(5): [SU(2) x C/"(l)]i x [SU(2) x U(l)] 2 , such that the so-called collective symmetry 
breaking prevents the Higgs from becoming massive when the couplings of one of the 
two gauge factors vanish. A discrete symmetry called T-parity [29] is then introduced, 
in order to reconcile the model with electroweak precision tests. It restores the custodial 
577(2) symmetry and, therefore, the compatibility with electroweak precision data is 
obtained already for quite small values of the NP scale, / > 500 GeV [61,62]. Another 
important consequence is that particle fields are T-even or T-odd under T-parity. The 
particles belonging to the T-even sector are the SM particles and a heavy top T + , while 
the T-odd sector consists of heavy gauge bosons W^, Zh, Ah, a scalar triplet $, an 
odd heavy top T_ and the so-called mirror fermions [63], i.e., fermions corresponding 
to the SM ones but with opposite T-parity and (9(1 TeV) mass. Mirror fermions are 
characterized by new flavour and CP-violating interactions with SM fermions and heavy 
gauge bosons, thus allowing significant effects in flavour observables [3-5,64-66] without 
new operators in addition to the SM ones. 

The formula for the CP-violating ratio e'/e of [23] in a generic model with new 
complex phases but no new operators, like the LHT model, generalizes as follows: 

e sm(/3 - p a ) 
with \ t = V*V td , (3 S = -1.3° and 

F £ ,(v) = P sin(/3 - p a ) + P E \E K \ sin /if 

+P x \X K \smp$ + Pr\Y K \sintff + P z \Z K \smpj[ , (3.2) 

where j3 is the angle in the unitarity triangle to be specified below (see Table 2). 

Pi are the same as in the SM while the short distance physics is now described by 
the loop functions 

X K = \X K \e i6 %, Y K =\Y K \e ie ?, Z K =\Z K \e^, E K =\E K \e^, 

(3.3) 

that are generalizations of the real valued SM loop functions X , Y , Z and E in (2.1) to 
the LHT model. Explicit expressions for X K) Y K and Z K have been obtained in [3]. The 



7 



function Ek can be found, in complete analogy to the functions Tp/ and T& governing 
the B — > X s 7 decay [64], by changing the argument of the SM E function and properly 
adjusting various overall factors. The result is given in Appendix A. The phases i3f 
entering (3.2) are then given by 

(3? = [3-p s -9? (i = X,Y,Z,E). (3.4) 

A comment on two approximations made above is in order. The first one concerns the 
contributions from the T-even sector to the functions Xk and Yk- In the calculation of 
these functions, the fermion mass on the flavour conserving side of the box diagrams has 
been set to zero, since in the case of semileptonic rare decays SM leptons are present. 
On the other hand, in the case of non-leptonic decays, such as Kl — > ttit, this mass 
cannot be generally neglected, as now up-type quarks, in particular the top quark and 
the heavy T + , contribute. However, it can straightforwardly be shown that including 
this difference results in the presence of a new operator [67] 

(sd) v _ A (bb)v_ A (3.5) 

at scales fx > nib, which is not contained in (2.1) and (3.1), (3.2). It is multiplied by the 
function 

S t = S (x t ) + Seven , (3.6) 

where So(xt) denotes the SM contribution and S even the heavy T + contribution. Below 
the scale fx = mb the b quark is integrated out, and therefore the operator in (3.5) 
contributes to e'/e only through mixing under renormalization. In the case of the SM, 
this contribution has been shown to be (9(1%) and therefore fully negligible [67]. As in 
the LHT model the dominant contribution to St comes from the SM part So(xt) [64,68], 
the accuracy of neglecting this contribution remains the same in the LHT model. 

The second approximation entering the above formula (3.2) concerns the T-odd sector 
and consists in neglecting the mass splittings of mirror quarks on the flavour conserving 
side of the box diagrams contributing to the Xk and Yk functions. We have checked, 
see also [3] , that the inclusion of these splittings affects the functions X K and Y K by at 
most 10%. As P x and P Y are much smaller than P and \Pz\, these functions do not 
play a dominant role in e'/e anyway and we can safely neglect also this effect in view of 
large non-perturbative uncertainties. 

In the LHT model, the first term in (3.2), which involves P and is dominated by the 
QCD penguin operator Q 6 , does not contain any NP contribution. On the other hand, 
the important negative last term involving P z and related to the EW penguin operator 
<2 8 can be strongly enhanced, when 9 Z 7^ 0, sin/3^ — 1 and \Z\ > Z (x t ). These 
conditions can indeed be satisfied, as found in [3] from a general scan over the three 
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generation mirror fermion masses and the six parameters (three angles 6f 2 , 0f 3 , 6^ and 
three phases Sf 2 , 8f 3 , S^) of the mixing matrix Vud- Thus, in this case, the suppression 
of e'/e through the enhanced electroweak penguin contribution must be compensated 
by the increase of the QCD penguin contribution Po or by decreasing the magnitude 
of the coefficient P z . This corresponds to the increase of P 6 and the decrease of R$, 
respectively. 

Clearly, as seen in the previous section, the result for e'/e is very sensitive to the 
actual values of the coefficients Pj. In the LHT model, in addition, there is a strong 
dependence on the phases fif. 

We conclude this section commenting on the origin of the correlations present in 
the LHT model between e'/e and rare kaon decays. They come from the simultaneous 
dependence of rare K decays and e'/e on the short- distance functions Xk, Yk and Zk- 
For instance, the branching ratio for K L — > tx q vv reads 

Br(K L -> 7i°uu) = k l ? 2 A a R 2 \X k \ 2 sin 2 p$ , (3.7) 

where [69] 

k l = (2.22 ±0.07) • 10" 

As, in the LHT model, there are also strong correlations between X K , Y K and Z K , 
in particular between their phases, it is evident that there will be a strong correlation 
between the CP- violating observables e'/e and Br(K L — > n°vv). 

The explicit expressions for Br(K + — > 7r + z/z/) and Bt{Kl — > 7r°£ + £~) in terms of Xk, 
Yk and Zk are more complicated than the one in (3.7). They are given in [3], to which 
we refer for details, and forecast as well correlations between e'/e and these decays. 

4 Numerical Analysis in the LHT Model 

In our numerical analysis in the LHT model presented below we have used for the 
determination of the CKM parameters, and in particular of lm(X t ), the tree level values 
of \V u b\, \V c b\, A and 7 given in Table 2, as the UTfit values obtained within the SM are 
clearly not valid in the LHT model. In obtaining the SM values of rare decay branching 
ratios in Table 3 below, however, we consistently used the determination of the CKM 
parameters within the SM. As a curiosity we remark that with the CKM values of 
Table 2, due to an increased value of Im(At) with respect to the UTfit determination, 
the SM branching ratios are higher than those given in Table 3 and their central values 



^-0.98, R t = ^P|-1.0. (3.8) 

Vcb \ V cd V cb\ 
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read 

Br(K L -> n°vvfiM = 4 -0 • 10~ n , fir(fsT + ->■ vr+z/z/)^ = 9.5 • KT 11 , (4.1) 

5r(fT L ->■ 7r°e + e - )^ = 3.8 • 10~ n , Br(K L ->• ttV + aOsm = 15 • 10_n ■ ( 4 - 2 ) 

However, such a procedure would not be fully consistent as the CKM values in Ta- 
ble 2 deviate significantly from the SM UTfit values: the reason is the so-called "sin 2/3 
problem" [70]. 



Gf 


= l.lboorf • 1U Lrev 


-2 


AM K = 3.483(6) • 10" 15 GeV 


M w 


= 80.425(38) GeV 




AM d = 0.508(4)/ps 


[72] 


a = 


1/127.9 




AM S = 17.77(12)/ps 


[73, 74] 


sin 


6 W = 0.23120(15) 


[31] 


S^ Ks = 0.675(26) 


[72] 


\v ub 


= 0.00409(25) 




F K \[B^ = 143(7) MeV 


[31,75] 


\Vcb\ 


= 0.0416(7) 


[72] 


F B d ^B Bd = 214(38) MeV 




A = 


\V US \ = 0.2258(14) 




Fb s \J Bb s = 262(35) MeV 


[75] 


7 = 


82(20)° 


[57] 


rn = 1.32(32) 


[76] 


m K { 


= 497.65(2) MeV 




% = 0.47(5) 


[77] 


m D o 


= 1.8645(4) GeV 




V2 = 0.57(1) 




m Bd 


= 5.2794(5) GeV 




rjB = 0.55(1) 


[78] 


™B S 


= 5.370(2) GeV 




m c = 1.30(5) GeV 




\ £ k\ 


= 2.284(14) • 10~ 3 


[31] 


m t = 161.7(20) GeV 





Table 2: Values of the experimental and theoretical quantities used as input parameters. 

The discussion of Sections 2 and 3 forecasts that in order to allow large enhancements 
of the rare decays Kl — > it !/!? and Kl — > tt°£ + £~ , the consistency with the data on e'/e 
requires i? 6 > R 8 . In Fig. 2 we show e'/e as a function of Br(K L — > n Q uu) in the 
LHT model for different values of (R e ,R$). To this end we have set A^jg = 340 MeV 
and performed a general scan over the parameters of the LHT model subject to present 
experimental constraints from K and B physics as discussed in detail in [3,64]. We 
compare the plot resulting from the general scan with the one obtained setting to zero 
two phases, 5f 2 and <5f 3 , of the Vud mixing matrix. 3 These two plots are significantly 
different, signaling that e'/e is quite sensitive to the new phases Sf 2 and #? 3 , whereas 
this sensitivity was much weaker in the case of rare decays discussed in [3]. This shows 
that e'/e is not only very sensitive to the values of the hadronic matrix elements but 
also to the new parameters of a given model. This fact could be used in the future to 

3 A detailed analysis of the number of phases in the mixing matrices in the LHT model has been 
presented in [71]. 
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Figure 2: Left: e'/e as a function of Bv(Kl — > tt°vi>) for different values of (Rq,R$) = 
(1.0,1.0) (red), (1.5,0.8) (green), (2.0,1.0) (blue). The shaded area represents the exper- 
imental result in (2.4) while the SM predictions are displayed by the black points. Right: 
Same as before, but with two phases (8f 2 o,nd 5 23 ) of the mixing matrix Vnd set to zero. 
Comparing the left and right plots, it is evident that e'/e turns out to be quite sensitive 
to these phases. 

efficiently exclude some portions of the parameter space provided the hadronic matrix 
elements will be brought under control. 

We observe that for (i? 6 , R$) = (1.0, 1.0) (red points), large enhancements of Br(Ki — > 
71"°^) over the SM value imply a strong suppression of e' / e relative to the data, and con- 
sequently in this case large enhancements of BriK^ — > 7r°^P) found in the LHT model 
in [3] are unlikely. The same applies to BriK^ — > tt°£ + £~). On the other hand, for 
(Rq, R$) = (1.5, 0.8) (green points) and (Rq, R$) = (2.0, 1.0) (blue points) the experimen- 
tal data for e'/e imply in the LHT model a significant enhancement of Br(K^ — > 7r°z/z/) 
with respect to the SM. 

As K~l — > n^vv and — > tt°£ + £~ are very strongly correlated with each other [3], 
also Br{Ki — > ti £ + £~) are predicted to be enhanced for (i? 6 ,i? 8 ) = (1.5,0.8) and 
(Rs,Rs) = (2.0,1.0). We summarize in Table 3 the three choices for (R 6 ,R 8 ) and the 
corresponding values of rare decay branching ratios that are compatible with the data 
for e'/e. 

In Fig. 3 we show the correlation between e'/e and K + — > 7i + uu that is significantly 
weaker than in the case of e'/e and K~l — > 7r°z/z/. In particular, we find that in the case 
(R 5 ,R 8 ) = (1.0,1.0), in which Br(K L -)■ -k°vij) and Br(K L -)■ n°£ + £~) are required 
to be SM-like, Br(K + — > 7t + uj?) can be largely enhanced relative to its SM value. A 
different behaviour is observed for the two other choices of (Rq,R$) considered by us. 
Here only enhancements of Br(K + — > 7i + uu) by at most a factor 3 are allowed. 

In order to understand better the pattern of enhancements of Br(K~L — > tx^vv) and 
Br(K + — > 7T + i>i?), we show in Fig. 4 the correlation between BriK^ — > 7i uu) and 
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Figure 3: Correlation between e'/e and K + — > tt + vv for different values of (R e ,R s ) = 
(1.0,1.0) (red), (1.5,0.8) (green), (2.0,1.0) (blue). The shaded areas represent the exper- 
imental results while the SM predictions are displayed by the black points. 




Figure 4: Left: Correlation between K + — > n + vv and Kl — > tx q vu without imposing the 
e' / e- constraint [3]. The shaded area represents the experimental result for Br(K + — > 
n + uu) [80] while the SM predictions are displayed by the black points. The Grossman- 
Nir bound [81] is displayed by the dotted line, while the solid line separates the two areas 
where BriK^ — > n^vv) is larger or smaller than Br(K + — > 7i + uu). Right: Same as 
before, but after imposing the constraint from e'/e with different values of (Rq,R$) = 
(1.0,1.0) (red), (1.5,0.8) (green), (2.0,1.0) (blue). 
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SM 


(1.0,1.0) 


(1.5,0.8) 


(2.0,1.0) 


Br(K L -> tt°vv) ■ 10 11 


2.7 ±0.4 


0.007... 9.5 


0.5... 43 


8.4. ..42 


Br(K+ ->■ tt+z/z/) • 10 10 


0.84 ±0.10 


0.09. ..5.7 


0.6... 2.3 


1.0... 1.8 


Br(K L ->■ 7r°e+e-) • 10 11 


o c 4 +0.62 
°-<- Kt -0.49 


2.7. ..4.7 


2.9. ..8.8 


4.2... 8.6 


Bt(Kl ->■ ttVV") ' 10 11 


1 41+0-28 


1.2... 1.8 


1.2. ..3.9 


1.8... 3.8 



Table 3: Choices for (R 6 , R 8 ) and the corresponding values of rare decay branching ratios 
that are compatible with the data for e'/e. The SM predictions [79] are also shown. For 
Br(K L — > ir £ + £~) we consider for simplicity only the case of constructive interference 
between direct and indirect CP-violation. 

Br(K + — > 7r + z/z/) in the LHT model without the e'/e constraint as obtained in [3], and af- 
ter the constraint from e'/e for different choices for (i? 6 , R 8 ) has been taken into account. 
We observe that setting (Rq,R 8 ) = (1.0, 1.0) basically selects the horizontal branch on 
which Br(KL — > 7r° z/z/) is SM-like but Br(K + — > can be strongly enhanced. The 

other two choices for (i? 6 ,i? 8 ) select the second branch on which Br(K^ — > tt uu) can 
be strongly enhanced but Br(K + — > 7i + uu) < 2.3 • 10~ 10 . 

5 Conclusions 

In this paper we have calculated e'/e for different values of the hadronic parameters 
(i?6, Rs) i n the LHT model and investigated the implications for rare decay branching 
ratios when taking the experimental data for e'/e into account. The main results of our 
paper are given in Figs. 2-4 and in Table 3 and can be summarized as follows: 

• For the values of hadronic parameters (i? 6 ,i? 8 ) — (1.0,1.0), for which (e'/e)sM 
agrees with the data, large enhancements of Br(K L — > ix^vv) and Br(K L — >■ 
ir°l + £~) relative to the SM are very unlikely. 

• On the other hand, for the values of hadronic parameters (Rq,R$) = (1.5,0.8) 
and (2.0, 1.0) chosen by us, the large NP contributions that are required to fit the 
experimental value for e'/e result in large enhancements of Br(Ki — > n vu) and 
Br(K L 7r°£ + £-) relative to the SM. 

• The correlation between e'/e and K + — > 7i + uu is much weaker and large departures 
of Br(K + — > n + vu) from the SM values are possible even for (i? 6 , R 8 ) ~ (1.0, 1.0), 
however, more modest enhancements are possible for the other choices of hadronic 
parameters, as seen in Figs. 3 and 4. 
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The main message of our paper is clear: without significant progress in the evaluation 
of Rq and R$ an d other less important hadronic parameters entering e'/e, the role of the 
data in (2.4) in constraining physics beyond the SM will remain limited. 
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A Explicit Formulae for the Function Ek 

In this appendix we give the explicit expression for the function Ek entering the calcu- 
lation of e'/e in the LHT model. The functions Xk, Yk and Zk have been calculated 
already in [3] in the context of rare K and B decays and can be found in that paper. 
The variables are defined as follows: 



m 



2 



- — S*. , z [ = a Zi with a = - — (« = 1, 2, 3) , (A.2) 

M W H tan vw 

At = V*V td , if ] = V™Vi? d (i = 1, 2, 3) , (A.3) 
and xl describes the mixing in the T-even top sector. 



E Q (x t ) 
Ek 



E 



odd 
K 



2, x?(15- lQxt + 4x?) , x t (18 - Ux t -x 



Q(i-x t y 



12(1 -x t f 



Eo(x t ) + E even + — E° K , 



E even = x\ — [E (x T ) - E (x t )] , 



4? 



(K) 



(A.4) 
(A.5) 
(A.6) 
(A.7) 
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